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ABSTRACT 



The Advanced Camera for Surveys (ACS) Virgo Cluster Survey is a large 
program to image 100 early-type Virgo galaxies using the F475W and F850LP 
bandpasses of the Wide Field Channel of the ACS instrument on the Hubble Space 
Telescope (HST). The scientific goals of this survey include an exploration of the 
three-dimensional structure of the Virgo Cluster and a critical examination of 
the usefulness of the globular cluster luminosity function as a distance indicator. 
Both of these issues require accurate distances for the full sample of 100 program 
galaxies. In this paper, we describe our data reduction procedures and examine 
the feasibility of accurate distance measurements using the method of surface 
brightness fluctuations (SBF) applied to the ACS Virgo Cluster Survey F850LP 
imaging. The ACS exhibits significant geometrical distortions due to its off-axis 
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location in the HST focal plane; correcting for these distortions by resampling the 
pixel values onto an undistorted frame results in pixel correlations that depend 
on the nature of the interpolation kernel used for the resampling. This poses 
a major challenge for the SBF technique, which normally assumes a flat power 
spectrum for the noise. We investigate a number of different interpolation kernels 
and show through an analysis of simulated galaxy images having realistic noise 
properties that it is possible, depending on the kernel, to measure SBF distances 
using distortion-corrected ACS images without introducing significant additional 
error from the resampling. We conclude by showing examples of real image power 
spectra from our survey. 

Subject headings: galaxies: distances and redshifts — galaxies: clusters: individ- 
ual (Virgo) — methods: data analysis — techniques: image processing 

1. Introduction 

The Virgo Cluster is an appealing target for the study of early-type galaxy formation 
and evolution because of its richness and its proximity. The installation of the Advanced 
Camera for Surveys (ACS; Ford et al. 1998) on the Hubble Space Telescope (HST) has 
provided the opportunity to observe Virgo galaxies with unprecedented spatial resolution 
in reasonable exposure times. The ACS Virgo Cluster Survey is a program to image 100 
early-type members of Virgo in the F475W (« SDSS g) and F850LP (« SDSS z) filters of the 
ACS Wide Field Channel (WFC). This survey aims to study the properties of the globular 
cluster systems associated with the program galaxies, examine the nuclear properties and the 
morphological structure of their central regions and explore the three-dimensional structure 
of the cluster using surface brightness fluctuations (SBF). A comprehensive description of 
the program is given in Cote et al. (2004; hereafter Paper I). In this paper, we discuss the 
feasibility of SBF measurements with the ACS. 

Over the last ten years, the SBF method has been successfully used to measure 
early-type galaxy distances out to ~ 7000 km/s from ground-based telescopes and HST 
(Luppino & Tonry 1993; Pahre & Mould 1994; Sodemann & Thomsen 1995, 1996; Ajhar 
et al. 1997, 2001; Thomsen et al. 1997; Tonry et al. 1997, 2001; Lauer et al. 1998; Jensen 
et al. 1998, 1999, 2001; Blakeslee et al. 1999b, 2001, 2002; Pahre et al. 1999; Mei et al. 2000, 
2001a,c, 2003; Liu et al. 2001, 2002; Mieske & Hilker 2003; Mieske et al. 2003; Jerjen et al. 
2001, 2003, 2004). The method was introduced by Tonry & Schneider (1988) and is based 
on the fact that the Poissonian distribution of unresolved stars in a galaxy produces fluctua- 
tions in each pixel of the galaxy image (see, e.g., the review by Blakeslee et al. 1999a). The 
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variance of these fluctuations is inversely proportional to the square of the galaxy distance 
and depends linearly on the mean flux of the galaxy To make these fluctuations constant 
on the galaxy image, the SBF amplitude is defined as this variance normalized to the mean 
flux of the galaxy in each pixel (Tonry & Schneider 1988). The absolute magnitude of the 
fluctuations depends on the age and metallicity of the stellar populations in the galaxy, 
which means that, in practice, the calibration of the SBF absolute magnitude for distance 
measurements depends on both the observational bandpass and the galaxy color. 

Different groups have quantified these dependencies for ground-based and HST 
photometric bands, empirically calibrating the dependence of the SBF magnitude on 
galaxy color (Tonry 1991; Tonry et al. 1997; Ajhar et al. 1997; Jensen et al. 1998, 2003; 
Ferrarese et al. 2000; Blakeslee et al. 2001; Liu et al. 2002). The dependence of SBF mag- 
nitude on color has also been studied extensively from a theoretical perspective using stellar 
population models (Worthey 1993; Buzzoni 1993; Liu et al. 2000; Blakeslee et al. 2001; 
Cantiello et al. 2003). At present, the absolute calibration is based on a comparison of SBF 
distances for galaxies with measured Cepheid distances. 

In using the ACS for SBF distance measurements, we are faced with two main ob- 
stacles. First, the ACS cameras suffer from significant geometrical distortions (Meurer et 
al. 2002). Correcting for these distortions involves resampling the pixels onto a rectified 
output grid which, depending on the choice of interpolation used for the resampling, may 
result in strong noise correlations that could bias the Fourier-space SBF measurements. It 
is therefore important to characterize this bias in order to ensure an accurate measurement 
of the Poissonian stellar population fluctuations. The second concern involves the proper 
calibration of the SBF method for the ACS passbands: a prerequisite for the measurement 
of SBF distances. Our survey will provide the first calibration of SBF distances in the ACS 
F850LP filter, which, as discussed in Paper I, is expected to be nearly optimal for the SBF 
technique. The first point will be addressed here, while the ACS SBF calibration will be the 
subject of a future paper. 

The plan of this paper is as follows. The following section discusses the ACS observations 
and the SBF data reduction procedures. In § 3, we investigate different interpolation kernels 
and demonstrate with simulations that our data reduction procedures allow us to make 
high-quality, essentially unbiased, SBF measurements. We conclude in § 4. 
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2. Surface Brightness Fluctuation measurements with the ACS 

A comprehensive description of the ACS Virgo Cluster Survey observations and initial 
data reduction can be found in Paper I and in Jordan et al. (2004; hereafter Paper II). We 
summarize here the main points. 

2.1. Observations and data reduction 

The ACS Virgo Cluster Survey consists of two-filter (F475W and F850LP) imaging of 
100 early-type galaxies in the Virgo Cluster, each one observed within a single HST orbit 
using the ACS/WFC. The data for each galaxy consists of five images: two exposure of 
375 sec in the F475W filter and three exposures in the F850LP filter (two 560 sec and one 
90 sec exposures, for a total exposure time of 1210 sec). 

The ACS detectors have strong geometrical distortions caused mainly by the off-axis 
location of the instrument (Meurer et al. 2002). This distortion causes the square pixels of 
the detectors to project to trapezoids with an area that varies across the field of view. By 
modeling and calibrating these distortions, it is possible to produce a geometrically correct 
image having a constant pixel area. As described in Paper II, we have used the drizzle 
software (Fruchter & Hook 2002) to combine the individual exposures of a given field in our 
program into a single geometrically corrected image for each filter. However, the interpolation 
necessary to generate this image can introduce correlations between output pixels. Since 
SBF is calculated as the variance of the Poissonian distribution of stars in each pixel, SBF 
measurements can be biased by this correlation. The amount of correlation that is introduced 
between nearby pixels depends sensitively on the choice of the drizzle interpolation kernel. 
If the images are not corrected for the distortion, and then not interpolated, the pixel area 
varies over the image (due to the distortion), also potentially biasing the measurement of 
SBF, since Poissonian fluctuations would be measured in pixels of different areas. Further, if 
there is dithering between multiple exposures, distortion correction is necessary in order to 
combine the images. We therefore measure SBF in corrected images, quantifying the bias in 
the SBF measurements by simulations. We describe the SBF reduction in the next section 
and discuss how we can obtain high-quality SBF measurements despite the geometrical 
distortion correction. 
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2.2. SBF measurements 

We measure the SBF signal in the F850LP data using the standard SBF extraction 
techniques (e.g., Tonry & Schneider 1988; Tonry et al. 1991). The SBF in the F475W data 
is too faint to measure accurately, but the galaxy color information from the two bandpasses 
will be an essential part of the ACS SBF calibration. The image processing that we perform 
in order to produce uniform, galaxy-subtracted images has been described in more detail 
by Paper II; here we summarize the main steps. We fit a smooth model for each galaxy 
in both filters using ELLIPROF (the isophotal fitting software that has been used for the 
SBF survey by Tonry et al. 1997) and then subtract the models from the original images. 
To eliminate large-scale residuals left from this subtraction, we run SExtractor and produce 
a smooth background model which we then subtract from the galaxy-subtracted image. 
Object detection is then performed on this residual image, using a detection threshold of 
five connected pixels at 1.5 a significance level (see Paper II for the details of these first steps 
of the data reduction). 

A catalog of sources is produced by matching the F475W and F850LP detections. This 
procedure is necessary to separate spurious detections (for instance, residual cosmic ray 
events or hot pixels) from real, faint objects. Red background galaxies, however, might 
produce a clear detection in F850LP only; we therefore added to the list of sources detected 
in both filters objects brighter than 23 mag in F850LP, but absent in the F475W frames. 
The source catalog was then trimmed by excluding objects fainter than a cut-off magnitude 
m cut . 

The final sources list thus obtained was used to create a "source mask," to remove the 
contribution to the image power spectrum from contaminating sources such as foreground 
stars, globular clusters, and background galaxies. We mask each source using an aperture 
with size equal to the maximum of 3 times the FWHM (Full Width Half Maximum) of the 
PSF (point spread function) (usually best adapted for compact objects) and the product of 
the Kron radius and the object major axis (best adapted to extended objects) as calculated 
by SExtractor. 

Detected sources with magnitude brighter than the cut-off magnitude m cu t are removed 
from the residual image by multiplying it by the source mask. The resulting image is then 
divided into concentric annuli centered on the galaxy center, and the power spectrum of 
each annulus is measured. To do this, we multiply the galaxy image by a total annular 
mask function defined as the product of the source mask function and the annular mask 
function. Assuming no noise correlations, the image power spectrum is the sum of two 
simple and distinct components: the flat, white noise power spectrum and the combined 
power spectrum of the SBF and the undetected faint sources, both convolved by the PSF in 
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the spatial domain. 

In the Fourier domain, the convolution by the PSF translates to a multiplication. We 
therefore fitted the azimuthally averaged image power spectrum P(k) to the function 

P{k) = P x E(k) + Pi , (1) 

where the "expectation power spectrum" E(k) is the convolution of the power spectra of the 
normalized PSF and of the mask function of the annular region being analyzed. In equation 1, 
Po represents the PSF-convolved component, while P\ is the white noise component which 
is modified by the distortion correction (see below). A high signal-to- noise composite z- 
band PSF was created using NGC 104 (47 Tuc) observations from HST calibration program 
GO-9018. After normalizing the total flux to 1 electron per second, its power spectrum was 
computed. While the PSF of the ACS WFC does show some spatial variation (Krist 2003), 
the variation in the z-band is not significant for our purposes as it is present on very small 
scales that do not significantly affect our power spectrum fits (see below). To verify this, 
we investigated the effects of using individual stars as PSF templates, and found that the 
typical variation in the fitted P was 2%. 

We wish to normalize the P component by the mean galaxy surface brightness in the 
region under consideration so that the SBF contribution to P is equal to the ratio of the 
second and first moments of the stellar luminosity function (e.g., Tonry & Schneider 1988). 
This is done by multiplying E(k), as defined above, by the square root of the galaxy surface 
brightness model for the annulus (Tonry et al. 1991). In practice, this is achieved by defining 
a "window function" as: 

W(pc) = M(x) x v/G(x) (2) 

where M(x) is the source mask and G(x) is the galaxy surface brightness model for the 
annulus. The 2-D expectation power spectrum is then 

£(k) = P PSF (k) ® iV(k) (3) 

where Ppsf(X) is the spectrum of the normalized PSF and -FV(k) is the power spectrum of 
the window function. The power spectrum is then azimuthally averaged, and a linear fit to 
equation (1) is made by minimizing the absolute deviation (Press et al. 1992) to derive P 
and P\. The uncertainty in the fit, including variations resulting from different choices for 
the adopted k range, is included in our final error estimates. A complete description of the 
errors will be given in the ACS SBF calibration paper. 

To our knowledge, all previous SBF measurements have been made without ever resam- 
pling the image pixel values (e.g., not correcting for distortions and allowing only integer 
pixel shifts); however, the instruments used for past SBF observations had significantly less 
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distortion than the ACS. The measurement of P and Pi in distortion-corrected images can 
be biased by the interpolation kernel used for the pixel resampling since the noise correla- 
tions (introduced by interpolation) change the shape of the power spectrum at the correlation 
scale. The choice of the interpolation kernel is therefore a critical aspect of the ACS data 
reduction, to which we now turn our attention. 



3. Simulations of power spectrum biases and results 

We have studied in detail the effects of four possible choices of the drizzle interpola- 
tion kernel on the SBF measurements: the square (Bilinear), gaussian (Gaussian), lanczosS 
(Lancsoz3; a damped sine function), and point (Nearest pixel) interpolation kernels. The 
sine function is the ideal low-pass filter, but it cannot be used by itself due to its asymptot- 
ically oscillating character. When used as an interpolation filter, it is usually multiplied by 
a window function. The 3-lobed Lancsoz-windowed sine function that we use is defined as 

sin(nx) sin(nx/3) i | g 

Lancsoz3(x) = { ~ x * x / 3 |X| (4) 

\x\ > 3 

To examine the effects of the interpolation kernel on the SBF measurements, we have 
generated two sets of simulated images. The first set contained simply a white noise compo- 
nent, while the second set consisted of model galaxies having total luminosities representative 
of the galaxies in the ACS Virgo Cluster Survey sample. The goal of these simulations is 
to study the effect of the pixel resampling under different interpolation kernels during driz- 
zling, rather than to simulate the ACS Virgo galaxy observations in full detail. Therefore, 
the simulated raw images do not have distorted pixels, but both sets of simulations are ge- 
ometrically transformed using drizzle in the same way that the actual distorted images are 
corrected in our ACS pipeline. 



3.1. Results from white noise simulations 

Sample images for the white noise simulations are shown in Fig. 1. It is apparent that 
both the Bilinear and Gaussian images show visible patterns of correlation. The Nearest 
and Lanczos3 kernels introduce less correlation. The Nearest interpolation leaves a large 
number of zero (uninterpolated) pixels in the output image and does not preserve relative 
astrometry at the pixel scale. This is because no real interpolation is done in this case; the 
input pixel value is simply assigned to the nearest output pixel. 
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To quantify the bias arising from the interpolation, we have calculated power spectra for 
the simulated images. For the white noise simulation, the original image has, as expected, 
a flat power spectrum. However, once the drizzling has been applied, both low and high 
wavenumbers are contaminated by the pixel correlations introduced by the interpolation 
kernels. We show in Fig. 2 our results for the four different kernels. In all cases, the original 
power spectrum is shown as a dashed line while the power spectrum of the drizzled image 
is shown as the continuous line. For comparison, Fig. 3 shows the difference between the 
power spectra measured in the original and drizzled images. 

As Fig. 2-3 reveal, the Bilinear and Gaussian interpolations introduce artificial pixel 
correlation on all scales, making it impossible to use these kernels for accurate SBF measure- 
ments. The Lanczos3 kernel depresses the power spectrum at high wavenumber (small pixel 
scales). The Nearest kernel does not bias in any significant way the image power spectrum, 
which remains flat (constant with k), even if the power spectrum is systematically lower 
than the original. From the noise simulations, both of the latter kernels can be used for 
SBF measurements in the k range in which the power spectrum remains flat ( i.e. provided 
the smallest pixel scales are not used when fitting the power spectrum with the Lanczos3 
interpolation). From the white noise simulations drizzled with the Lanczos3 kernel, we have 
measured the pixel scale range in which the power spectrum is flat to within 3 a. This 
is true for scales between 3 and 20 pixels, corresponding to fractional wavenumbers in the 
range 0.05 < k < 0.33. On smaller pixel scales, the original white noise power spectrum is 
significantly depressed by the noise correlation. 



3.2. Results from the galaxy simulations 

The galaxy simulations have been performed following the method described in Mei et 
al. (2001b). With the galaxy magnitude and effective radius specified, the surface bright- 
ness at each pixel is calculated using a de Vaucouleurs (1948) profile and Bruzual & Chariot 
(2003) stellar population models corresponding to solar metallicity and an age of 12 Gyr. 
The results are not sensitive to the use of the de Vaucouleurs profile, since the SBF mea- 
surements are done on the galaxy model-subtracted images. The isophotal model used for 
the galaxy subtraction does not assume the de Vaucouleurs form (see §2.2). Just as in real 
data, mismatches between the simulated galaxy images and isophotal models will produce 
large-scale spatial signatures in the power spectra that do not affect our fitting scales. Simi- 
larly, the precise age and metallicity of the stellar population model is not important, except 
insofar as they precisely specify the SBF magnitude for the simulations. 

To each pixel in the simulated images, we add a Poisson variance calculated from the 
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number of stars contributing to the pixel, Poissonian noise from the galaxy and the sky 
(22 magarcesc -2 ), and a read noise of 7 e~. The images are then convolved with the PSF 
and drizzled using three different interpolation kernels: Gaussian, Lancsoz3 and Nearest. 
Since the results from the white noise simulations show that the Gaussian and Bilinear 
kernels introduce similar strong pixel correlations, we use only the Gaussian to quantify 
this effect. Moreover, in the larger context of our survey, this kernel allows a more faithful 
reconstruction of the galaxy centers. Galaxy simulation images are shown in Fig. 4 and 
Fig. 5. The left panels show the simulated galaxies and the right panels show the residuals 
following subtraction of smooth isophotal models. As for the white noise simulations, the 
Gaussian interpolation shows visible patterns in all the residuals. The Nearest interpolation 
shows visible residual patterns in the centers of the galaxy simulations where the intensity 
profile is steepest, especially for the brighter galaxy. 

To quantify the effect of the interpolation on the SBF measurements, we have simulated 
100 galaxies for an assumed Virgo distance of 16 Mpc, and for three galaxy magnitudes 
corresponding to: 1) the brightest galaxy in our sample (Mb = —22 mag); 2) the faintest 
galaxy (Mb = —15 mag); and 3) an intermediate galaxy (Mb = —17 mag). The B — V 
(Vega) color from Bruzual & Chariot (2003) models is 0.97 mag, and the V — F850LP 
(AB) color is 0.9 mag. Each of these simulations has been drizzled with the three different 
kernels. We then applied the SBF reduction procedure, as described above, to the simulated 
galaxy images to obtain SBF magnitudes, and then distances from the input SBF absolute 
magnitude (M fs50lp = —1.42 from the Bruzual & Chariot stellar population model). 

Histograms of the resulting SBF distances for individual annuli of the 100 galaxy image 
realizations are shown in Fig. 6 and summarized in Table I. SBF magnitudes and distances 
are shown when the fit is performed in two different ranges: (1) between 3 to 20 pixels, and 
(2) between to 20 pixels, i.e., out to the maximum wavenumber. The fit was not extended 
beyond 20 pixels because, in the real data, larger scales will be biased by the subtraction of 
the large-scale residuals (see § 2.2). 

The Gaussian filter introduces a bias in SBF distance measurements of 6 to 20%, com- 
pared to a bias of 0.1 to 4% for the Lancsoz3 kernel. The Nearest filter tends to underestimate 
the distance in the innermost annuli of bright galaxies. As seen in Fig. 6, when the Nearest 
filter is used, the distance histogram of the brightest galaxy shows a bimodal distribution, 
with peaks at ~ 11 Mpc and « 15 Mpc, corresponding to the innermost and outermost 
annuli, respectively. This causes a bias of 20% on the total SBF distance of bright galaxies. 
For fainter galaxies, the biases due to this filter are between 1 and 2 %. However, contrary 
to expectations from the white noise simulations, when the Lanczos3 filter is used, galaxy 
distances do not significantly improve if we use a maximum wavenumber corresponding to 3 
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pixels. This is because of the reduced leverage for determining the Pi value in equation (1). 

We have chosen to use the Lanczos3 kernel for our final SBF analysis to insure a homoge- 
neous data treatment between faint and bright galaxies, and because it preserves astrometric 
fidelity. This has the further advantage that the same images can be used for fitting globular 
cluster structural parameters. However, since the Lanczos3 kernel is not well adapted to 
isophotal analysis (e.g., negative pixel values occur near saturated or repaired pixels), we 
have produced a second set of images with a Gaussian kernel, used in the isophotal analysis 
of the program galaxies. As a visual comparison of the effects on real data, the top panel of 
Fig. 7 shows a corner of the M87 image drizzled with the Gaussian (left) and the Lanczos3 
(right) kernel. The lower panel compares the respective weight images for the two kernels as 
derived from the drizzle interpolation (see Paper II). 

Fig. 8 shows simulated galaxy power spectra in an annulus with inner and outer radius 
of, respectively, ~ 1" and ps 6". On the left, the spectrum of a galaxy with an absolute B- 
band magnitude of —22 mag is shown before (top) and after (bottom) drizzling; on the right, 
the same comparison is given for a galaxy with absolute 5-band magnitude of —17 mag. 
We have fitted the power spectrum of the original image over the entire scale range and the 
power spectrum of the drizzled image for scales between 3 and 20 pixel. For this particular 
choice of annulus, the distance values between the original and drizzled simulated images 
agree to 4%, with a corresponding difference in SBF magnitudes of ~0.05 mag. In Fig 9 we 
show NGC 4472 power spectrum for three innermost annuli with respective inner and outer 
radius of 1"6 and 372, 372 and 674, and 6"4 and 976. 

4. Conclusions 

We have described the SBF data reduction procedures for the ACS Virgo Cluster Survey, 
paying close attention to the feasibility of SBF measurements with the ACS/WFC instru- 
ment. Although the ACS images exhibit strong geometrical distortions that could bias SBF 
measurements, we have shown with the aid of simulated galaxy images having realistic noise 
properties that it is possible to measure accurate SBF distances by correcting for biases, 
taking care to mask contaminants (i.e., globular clusters, foreground stars and background 
galaxies), and using an appropriate (Lanczos3) interpolation kernel. Future papers in this 
series will derive a new z-band SBF calibration for our sample, present the final SBF dis- 
tances and their uncertainties, compare the SBF and globular cluster luminosity function 
distances, and examine the three-dimensional structure of the Virgo cluster. 
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Table 1: Measurements of SBF Magnitudes and Distances from Drizzled Simulated Images. 
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(mag) 
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-22 


Original 


0-20 


29.57 


0.08 


15.9 


0.6 


0.6 




Gauss 


0-20 


29.46 


0.08 


15.0 


0.6 


6 




Lanczos 


0-20 


29.60 


0.08 


16.0 


0.6 


0.1 




Nearest 


0-20 
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0.66 
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24 
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Original 
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0.08 


16.1 
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Fig. 1. — Results from "drizzling" white noise simulations with different interpolation ker- 
nels. On the very top at the center is shown the white noise simulation before drizzling. The 
central two images have been drizzled with the Bilinear (left) and Gaussian (right) kernels. 
The lower images have beens drizzled with the Lanczos3 (left) and Nearest (right) kernels. 
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Fig. 2. — Noise power spectra when different interpolation kernels are used in the image 
drizzling. From top to bottom, the four panels refer to the Bilinear, Gaussian, Lanczos3, 
and Nearest kernels. The wavenumber k corresponds to the inverse of the pixel scale. The 
continuous lines are the spectra after drizzling, the dashed line is the power spectra of white 
noise before drizzling. 
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Fig. 3. — Differences between the original and drizzled noise power spectra are shown as 
continuous lines for the same four cases as in Fig. 2. From top to bottom, the panels refer 
to the Bilinear, Gaussian, Lanczos3, and Nearest kernels. The dashed line is the case where 
the difference is zero. 
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Fig. 4. — Results from galaxy simulations with different interpolation kernels. On the very 
top is shown a simulation (before drizzling) for a galaxy (ss 30 " x 20 ") with Mb = —22 
(left), and the residual after subtracting a smooth galaxy model (right). The residuals show 
fluctuations from Poissonian noise, read out noise and the PSF-convolved Poissonian stellar 
fluctuations. In the center images, the same for the same galaxy image drizzled with a 
Gaussian and a Lancoz3 filter. In the bottom, the same, drizzled with a Nearest filter. 
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Fig. 6. — The histograms of the measured SBF distances are shown for the three different 
magnitude simulations. For each of the 100 simulations at each magnitude, SBF magnitudes 
and distances have been measured in four different annuli having inner and outer radii of l'/6 
and 3"2, 3"2 and 4"8, 4"8 and 6"4, and 6"4 and 8". From left to right we show results for the 
original galaxy simulations, and, respectively, for the Gaussian, Lanczos3, and Point/Nearest 
drizzled images. From top to bottom the galaxy absolute B magnitude is —22 mag, —17 mag 
and —15 mag. The solid histograms correspond to results obtained with a fit on the entire 
scale range, the dashed ones to fits obtained in the range between 3 to 20 pixels. The dashed 
vertical line corresponds to the input distance of 16 Mpc. 
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Fig. 7. — Top: a 0!5 x 0f7 section of the M87 image is shown when the Gaussian (left) 
and Lanczos3 (right) drizzle interpolation kernels are used. Bottom: the full drizzle weight 
images obtained for these interpolations (including weight variations from rejection of pixels 
affected by cosmic rays and detector defects), again with the Gaussian results on the left 
and the Lanczos3 results are on the right. 
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Fig. 8. — Original (top) and drizzled (bottom) power spectra for the galaxy simulations 
when a Lanczos3 kernel is used. The results are for an annulus with inner and outer radius 
of respectively ~1 and ~6 ". On the left, the spectrum of a galaxy with Mb = —22 mag 
is shown before (top) and after (bottom) drizzling; on the right, the same thing is shown 
for a galaxy with absolute Mb = —17 mag. The diamonds are the image power spectra, 
the continuous line the fit of the power spectrum P(k), the dotted line is Po x E(k) and the 
dot-dashed line indicates the fitted P\. 
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Fig. 9. — VCC1226 (M49, NGC 4472) power spectra for the three innermost annuli. From 
top to bottom, the power spectra were calculated using annuli having inner and outer radii 
of 1'.'6 and 3'.'2, 3'.'2 and 6'/4, and 6'.'4 and 9'.'6. The diamonds are the image power spectrum, 
the continuous line the fit of the power spectrum P(k), the dotted line is P x E(k), and the 
dashed line is P±. 



